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We study analytically the distribution of fluctuations of the quantities whose average yield the 
usual two-point correlation and linear response functions in three unfrustrated models: the random 
walk, the d dimensional scalar field and the 2d XY model. In particular we consider the time 
dependence of ratios between composite operators formed with these fluctuating quantities which 
generalize the largely studied fluctuation-dissipation ratio, allowing us to discuss the relevance of the 
effective temperature notion beyond linear order. The behavior of fluctuations in the aforementioned 
solvable cases is compared to numerical simulations of the 2d clock model with p = 6, 12 states. 
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I. INTRODUCTION 



While in equilibrium systems the correct weighting of fluctuations has disclosed the basis of modern statistical 
mechanics, uncovering the properties of non-equilibrium fluctuations remains one of the most challenging and far 
reaching open questions of statistical physics. However, in spite of recent important developments explicit calcu- 
lations, especially for interacting systems, are limited to a few cases Q. While for non-equilibrium stationary states a 
major advance has been the recognition of a general symmetry of the probability distribution of certain observables, 
described by the so-called fluctuation theorems Q, the behavior of fluctuations in slowly relaxing systems such as 
those undergoing critical dynamics, phase-ordering kinetics, glassy evolution, forced relaxation and others, is much 
less understood. Traditionally these systems have been characterized in terms of the space and time dependence of 
two-point correlation functions and responses (see, e.g., the review articles 0-0] )■ It has been only relatively recently 
that interest has moved to the study of fluctuations of these quantities and, in particular, their characterization via 
higher order correlation and response functions 043- These properties should, obviously, give us a more detailed 
description of the behavior of these highly non-trivial processes. A large activity around the analysis of fluctuations 
in glass forming liquids and other similar materials exists [l3[ . 

In this paper we study analytically the distribution of fluctuations of some quantities the average of which yield 
the usual two-point correlation function and linear response in three unfrustratcd models: the random walk, the d 
dimensional scalar field and the 2d XY model. In its structure this paper follows the presentation in [l4| where 
aging and the usual fluctuation-dissipation (FD) relation were studied in these same three models with the purpose of 
highlighting the fact that neither disorder nor glassiness are needed to obtain non-trivial dynamic features. Composite 
operators formed with the above-mentioned fluctuating quantities are related to higher order correlation and response 
functions 0, [ll|, EH Ell , and ratios between such composite operators can then be defined which generalize the largely 
studied FD ratios. The analysis of these ratios in the asym pto tic time limit inform us, in particular, on the relevance 
of the effective temperature notion 16] (for reviews see [17 19]) beyond the linear order. These analytical studies are 
complemented by numerical simulations of the 2c? clock model with p — 6 and p — 12 states, where a similar behavior 
is found. 

The paper is organized as follows. In Scc.|TT]wc define the quantities that we will compute. Section Hill is devoted to 
the analysis of Brownian motion in the over-damped limit. In Sec. IIVI we solve the Langcvin dynamics of the scalar 
field model in d dimensions. In Sec. [V] we study the evolution of the clock model and its limit, the XY model, in two 
dimensions where a Kosterlitz-Thouless (KT) phase exists in which the models are critical. Finally, in Sec. IVII we 
present our conclusions. 



II. FLUCTUATING QUANTITIES 



We wish to characterize the statistics of fluctuating quantities the average of which yield the two-time correlation 
and linear response. The choice of these fluctuating quantities is not unique and their meaningful definition depends 
on the problem of interest, as we will discuss below. In the simpler cases (as the random walk of Sec. IHII or the scalar 
field considered in Sec. IIV[) one is usually interested in the correlation and response of a field whose dynamics is ruled 
by a Langevin equation of the type 

T being an effective Ginzburg-Landau free-energy and £ a thermal noise with the usual properties (£(x,t)) = and 
(£(x, t)£(x' , t')) = 2TS(x — x')8(t — t') (we set Ub = 1 and we absorbed the coefficient T with a redefinition of time). 
This equation is complemented by the choice of an initial condition 4>{x,t = 0). 

In this stochastic dynamic equation there are two sources of fluctuations: the initial condition, if taken from a 
probability distribution, and the random force mimicking the thermal noise. In the process of computing a response 
function, the applied perturbation may be an additional source of fluctuations, if it is taken to be random. In the 
analytical approaches contained in this paper we will choose to work with fixed initial conditions and we will therefore 
freeze the first source of fluctuations mentioned above. The effect of random initial conditions will be considered in 
Sec. lVB~2l for a discrete model, the clock model, which cannot be described in terms of a Langevin equation and must 
be studied numerically. 

Let us now introduce the two-time average quantities we will be interested in, the correlation and the linear response 
function, and their fluctuating parts. The two-time correlation function is defined as 



C = C(x, x 1 ; t, t') = ((f>(x, t)4>(x', t')} , 



(2) 
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where x and x' are two generic points in space and t and t' two generic times. From this expression the definition of 
the fluctuating part of the correlation as 

C = C{x,x';t,t') = 4>{x,t)4>(x" ,t') , (3) 

is meaningful, since taking the average of the quantity in Eq. ([3]) one immediately recognizes the (averaged) correlation 
C = C(x, x'; t, t'). Although the definition (J3)) is not unique, since other different fluctuating quantities may share the 
same average, for the correlation function this appears as the most natural choice. Let us notice that not only the 
average of 4>(x, t)<fi(x' , t') yields C(x, x'; t, t'), but composite operators (or moments) formed with this fluctuating field 
provide higher order correlation functions. This observation, which may seem quite trivial for the correlation function, 
is not so obvious when dealing with the response function, for which the definition of a fluctuating part is less straight- 
forward. The usual perturbation used to compute it is a field, say h, possibly drawn from a probability distribution, 
that couples linearly to the variable of interest, say cf> itself, in such a way that T — > T — J d d x h(x,t)<f>(x^,tY The 
perturbed Langevin equation then acquires an additional additive force. From here one can easily prove pL [Til fl9|. 
for any J 7 , that the averaged linear response function is simply related to the correlation between <f> and the noise £: 



R(S,x';t,t')- 6{<i,h ^ t)) 



Sh(x',t') 



= ^{mms\t')). (4) 

h— 



This relation gives us two fluctuating fields, |^rpyU=o an( l <t>(x> ,t')/(2T), whose noise averages yield the 
linear response function, as natural candidates to define the fluctuating part. The question then arises as to which is, 
between these two, the most interesting quantity to consider. In fact, although these objects have the same thermal 
average, they may have different fluctuation spectra and higher-order correlations. In [l2| we showed that the two-time 
fields appearing in the Martin-Siggia-Rose-Jenssen-dcDominicis dynamic generating functional are naturally related 
to <j)(x,t)^(x',t')/(2T). This suggests that the quantity ()>(x,t)£(x',t')/(2T) may play a special physical role and is 
worth of being studied. Indeed, special and interesting properties are exhibited by this fluctuating quantity in the 
case of aging spin-glasses (l2| . Similar considerations have led to use this same fluctuating quantity in ferromagnetic 
systems in [JJJ. Furthermore, it must be recalled that composite operators formed with the field cf>(x, t)^(x' ,t')/(2T) 
are directly related to higher order response functions 0, [HI [TH, Elf , a property analogous to the one discussed above 
for the fluctuating part of the correlation function. This allows one to define ratios between composite operators of 
the fluctuating parts of response and correlation (see Sec. Ill B[) . which are natural generalizations of the usual FD 
ratio and, by studying their behavior, to discuss the relevance of the notion of an effective temperature beyond linear 
order. Finally, restricting to the quadratic models considered analytically in this paper, the quantity gf^ff)^ \ h=o 1S 
deterministic, while cj)(x,t)£(x' ,t')/ (2T) has a non-trivial fluctuating pattern. Indeed, considering as an example a 
free-energy functional of the form J- = J d d x (V</)) 2 /2 (but the argument is general for any quadratic form), one has 
that the solution of Eq. (TT]) is 

(j> h (k, t) = 0(fc, 0)e" fc2t + / dt' e-^*-*'' [£(fc, t') + h{k, t')} (5) 
from which one immediately concludes that 



84>h{k,t) 
5h(k',t>) 



is a deterministic function, independent of the initial condition, the applied field and the thermal noise realization. 
Instead, (f>(k,t)£(k',t')/(2T) fluctuates and its higher order moments are non-trivial, as we will show in the following. 
In conclusion, we define the fluctuating part of the linear response as 

R = R(x,x';t,t') = ^x,t)^^- , (7) 

which, together with Eq. ([3]), provide the definitions of fluctuations of the two-time quantities that will be adopted 
throughout this paper. 

In the following we will also consider the time-derivative of C: 

~ fir 1 

d t ,C=^ F (x,x , ;t,t'), (8) 
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and the integrated fluctuating linear response function 

ti = $(x,x';t,t'X) J* %f^ dz - (9) 

With the choice t" — t the average of the latter quantity provides the dynamic (sometimes denoted as zero field 
cooled) susceptibility 

X (x, x"; t, t 1 ) = {x{x, x 1 ; t, t')) = j ds R{x, x'; t, s) . (10) 
For the 2d XY model we will define the relevant fluctuating quantities in Sec. [V] 

A. Joint probability distribution 

In simple models with a quadratic Hamiltonian H, the probability distribution of the fluctuating quantities intro- 
duced above can be explicitly exhibited, and a full characterization of the fluctuation spectra can be given (for the 2d 
XY model an equivalent delineation will be provided by studying the moments of any order in Sec. [Vj)- Let us start 
with the simplest case of a single time-dependent function <f>(t) (i.e. we consider a zero space dimensional problem; 
namely, there is no x dependence). It is convenient to introduce the quantities 

§{t",t') = <t>(t")-4>{t'), (ii) 

and 



t" 



dz . (12) 



In terms of these quantities the variation of the correlation AC = C(t,t") — C(t,t') and the integrated response 
function fl read 

AC(t,t",t') = 4>(t)W,t'), (13) 

and 

^. .. , N 6(t)E(t",t') , , 

£(MV) = - ' 2 ^ - ■ (14) 

(With an abuse of language we omit to call "fluctuating" the hat quantities.) Notice also that the differential quantities 
dfC and R can be obtained from these forms as dfC = lim,5_j.o AC/ 5 and R = lim^o /V<5, where S = t" — t' . The 
joint probability distribution V(<f>, $, S) of having </> at time t, 4> at times t' and t", and 2 also at times t' and t" is 
Gaussian and reads 



$ H ) G- 1 $ 



(15) 



iAC i/i 



P((/>,$,S) = (2tt)-' |G|-5exp 
where G is the matrix of correlations: 

G= \\AC ($ 2 ) \m I , (16) 

\G\ its determinant, {(j) 2 ) = C(t,t) is the field correlator, ($ 2 ) = C(t",t") + C(t',t') - 2C(t",t') is the field "displace- 
ment", AC = (AC(t, t", t')) = C(t,t") — C(t,t') and /i = (p,) = J dz R(t,z) in the time-integrated linear response. 
For white noise (S 2 ) = 2T(t" — t'). Let us observe that, by choosing </>, S as arguments of V one has the advantage 
of having all finite entries in G [this is not ensured if one uses, for instance, <fi(t), 4>(t') and £(£')]. 
In terms of V, the joint correlation-response probability distribution P(AC,/t) reads 

P(AC,-p)= J dcj)d<PdZT((j),<P,Z)5[<f)<P-eAd}5[(j)Z-2T6'fi} . (17) 
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Notice that we have introduced the parameter e which allows us to use a single form for the probability P(AC, ju) of 
the integrated quantities, by simply setting e = 1 in Eq. ([IT]) , and the probability P(dfC, R) of the differential ones 
by letting e = S and taking the limit S — > 0, namely 



P(d t >C,R) = 



lim P(AC,/2) 



5^0 



AC=d t , C 



(18) 



Using the integral representation S(x) = (2tt) 1 dr\ e ir,x of the Dirac function one arrives at 

roo+iv* drj 



P(AC,j£) = (2tt)-5 |(?|-3 d<j)d<I>dZ 



x exp 



oo+iA* j\ ^ ^ 

" A e - I (i)eAC+2ATe ( i) 



27r — oo+iA* 27r 



$ S)B. 



-l 

'/A 



(19) 



with 



2 





' 2 



(20) 







Notice that, following standard techniques [2(|, the integration paths in Eq. (|19p have been deformed by the arbitrary 
shifts rj* and A*. This docs not change the value of the integral since it merely amounts to introducing an additive 
term 77* (</><!> — eAC) + A*(</>S — 2Te/2), which vanishes due to the (^-constraints in Eq. (fl"7|) . in the argument of the 
exponential on the right-hand-side (rhs) of Eq. (|19[) . In turn, a properly deformed integration path may render the 
integrations over d(f)d$>d£ convergent in the cases in which (like in Sec. I IV A[) the result over the orginal route is 
divergent (when the integrations over d(j>d$>dE are taken before those over drjdX). 
Performing the Gaussian integrations in Eq. (|19|) one obtains 



P(AC,£) =|G|- 



00+171 dq 

— 00+171* 



oo+ir, ^ 

00+2A* 



oo+iX* j\ 

21 e -i(r,eAC+2AT eM ) \B x \? 



— oo-\-i\* 



2^ 



-i(r) eAC+2\Te p.) 



1 - -irjAC - -iA(27» 



■p\[2T^C 2<0 2 )(<I>~>] - ^ ?? 2 [(AC) 2 - 4(0 2 ><$ 2 ) 



1 

To 



A 2 [(21» 2 -4(,/> 2 )<S 2 >] 



(21) 



Note that the three time-dependencies (t > t" > t') enter only via averages. We will explicitly compute P in a specific 
zero-dimensional model, the random walk, in Sec. MI Al 

Let us now generalize what done insofar to a scalar field 4>(x, t) defined onad+1 dimensional space. Taking into 
account the space dependence the definitions (|11|) and (|T2|) arc replaced by 



$(x',t",t') = <j>{£ ,t") - <t>{af J) , 



and 



Z(x',t",t' 

The fluctuating quantities we are interested in are 

AC(r = x- x',t,t",t') = 



t" 



£,(x", z) dz 



dx cj>(x,t)<S>(x" \t" \t') 
dk 



(22) 



(23) 



(2tt) c 



(24) 
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and 



$(r = x-x',t,t",t') = 2y I dxcl>(x,t)E{x',t",t') 



(25) 



where <p(k,t) is the k component of the Fourier transform [5a | of the field <j)(x,t), and similarly for <3> and H. From 
here onwards r = x — x' . When the k components are independent (i.e., for a quadratic Hamiltonian) , the joint 
probability P(</>, <3>, S) is factorized as 



V(cf>^,E) = H(2nyi \G\-i exp 



t 



(26) 



where ( $ E ) and $ are evaluated at the wavevectors fc and —k, respectively. G is the matrix of the 



/c-component correlations 



\AC ±m 



2" /l ~2 C<J>H 

with AC = C(k,t,t") - C(k,t,t'), where 

(2ir) d 5(k + k')C(k,t,t') = (<j>(k,t)<p(k',t')), 



(27) 



(28) 



is the usual two-time structure factor. The other elements of the G matrix, which as in the zero-dimensional case 
are all finite, are the correlators (27r) d £(fc + k')]l(k, t, t",t') = (<j)(k, t)E(P, t", t')), (27r) d <5(fc + k')C<s>s(k, t, t", t') = 
($>(k,t" ,t')E(k' ,t" ,t')), and similarly for C^, C$$, and Chh. 

Starting from this, the joint probability distribution P(AC, fl) of the real-space quantities defined in Eqs. (|24[) and 
([25]) can be straightforwardly obtained proceeding along the same lines as for the zero-dimensional case. We first 
introduce the matrix B^x which is related to G by 

/ 



B- x 1 (k,f=x-x') = G- 1 (k) + 



Q x') _^ e ik(x-x') 



\ 



— in. e ik(x-x') 

:\ e ik(3—x') 
t 2 e 



(29) 



and is the generalization of Eq. (|20|) . We next obtain 



P(AC, V) 



— oo+ir/* 27T J—oo+iX* ^TT 



exp 



, 2^ 



°°+ lA * dX 

—— exp 



oo + iA 



, 2tt 



xcxp -| 



dk 



-i(r]eAC + 2ATe£) 



-i(r]eAC + 2XTefl) 



^Inll-^e^AC 



V 
~2 



dk 
(2n) d 



In 



I^A(fc,r)| 
\G(k)\ 



(2tt) c 



\iXe liF {2TJi) - lriXe 2i ^[2TJlAC - 2C H C* B ] - 
z o 

l^e^KAC) 2 - 4GmC**] - -^X 2 e 2 ^[(2Tji) 2 - AC^C SS ] 



(30) 

where V is the volume of the system, which extends Eq. (f2~Tj) to the finite dimensional case. This form is totally general 
(for a quadratic Hamiltonian). Note that the three time dependencies (t > t" > t') still enter only via averages, as in 
the zero-dimensional case, while the r- (actually r-) dependence is explicit. We will discuss the properties of P in a 
specific d-dimcnsional model, the scalar field, in Sec. IIV Al 
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B. Composite operators 

We now go on by using the generic notation and we build averages of composite fields of the form 

n m 

V (n,m) = { g t , d . JJ £ . } (31) 

i— m+1 i = l 

where we have used the shorthand 

C'i = C^Xi^x^tijtj,) (32) 

and similarly for Ri, where a% and x\ [ti and i£) are two generic space positions (times), with U > t\. 

It is also interesting to use the quantity in which the fluctuating time- variation of the correlation and linear response 
have been time-integrated: 

n m 

c (n ' m) =( n d ' n>> ( 33 ) 

i—m+l j — 1 

Moments are obtained from these quantities by subtracting a suitable disconnected part, as will be discussed in 
Sec. IV ATI The averaged correlation, linear response and the integrated linear response are simply C^ 1,0 \ XK 1 '!) ancl 
£(M) ^ respectively. Analogously, for n > 1 the quantities C^ n '°' > (or D^™' )) are higher order correlation functions 
(or their time derivatives) and, similarly, as discussed in [Til El 01, c(n ' n) and v(n ' n) 

are related to higher order 

response functions (time integrated or impulsive, respectively) [561 ]. 

Notice that the quantities in Eqs. pip and Q33p cannot be obtained from the joint probability distribution computed 
in Sec. Ill Al indeed, in general, they involve correlations and responses evaluated at different space-time variables 
Xi, x\\ U, t\ for any i, whereas both AC and /2 in Eqs. (|24p and (|25p are considered with the same space-time arguments 
(namely, x,x' and t,t',t"). 

Let us now define the generalized FD ratio as 

2?( n ' m ') ' 

Letting n = 1, m — 1, m' = one recovers the usual FD ratio X. These quantities do not put higher order FDTs to 
the test directly as these are complicated functions, involving several terms [si, [Hj], but they do evaluate whether all 
these terms scale in the same way thus allowing for the existence of an effective temperature taking finite values over 
non-trivial time regimes. 



III. BROWNIAN MOTION 



The ovcrdamped Langevin dynamics of a Brownian particle is ruled by Eq. ([T]) with a single function <f> (i.e. there 
is no dependence on x) and % = 0. In this case 4> should be interpreted as the position of a Brownian particle on a 
line. The extension to the case of a vector (j) with d components, describing diffusion in d dimensions, is trivial and 
the basic results remain unaltered. 



A. Joint probability distribution 

In this simple problem one trivially has 4> = 5. The joint probability distribution ([21]) of AC(i, t"t' ) and Jl(t, t" , t') 
then reduces to 



p(Ad,$) = 




x , 1 = (35) 

^1 - 1% + A)AC - ^(77 + A)2[(AC)2 - 4(^)($2)] 

with AC = 2T(t" - t'), (</) 2 ) = 2Tt and ($ 2 ) = (S 2 ) = 2T(t" - t'), where we assumed that t > t" > t' and set 
if = A* = since every integral is convergent. 
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This form is symmetric under the exchange AC 2TJi, thus indicating that these two quantities are equally 
distributed. Indeed, since $ = 5 one has AC = 2T/2 and, therefore, the fluctuations of the composite fields whose 
averages are the correlation and linear response are just identical in this case. As a consequence, all the composite 
operators of Eq. (|3"Tj) with the same value of n [with (ij,fj) = (t,t') Vi] scale in the same way (that is to say, they 
have the same dependence on the times t, t'). Moreover, introducing uj = r/ + A one can explicitly integrate expression 
(1351) and find 



P(AC ,jZ) = 6[ e(2T/i - AC) 



-iueAC 



2n 



y/l - \iuAC - > 2 [(AC) 2 - 4<0 2 }<$ 2 >] 



In the above integral, the integrand has two branch points at oj = u± , with 

±2 v /(0 2 )($ 2 )-AC 



lj± = 4 ■ 



(AC) 2 - 4(0 2 }($ 2 



(36) 



(37) 



Performing the integral we obtain 



4ACeAC 



P(AC, m) 



4e 



4(0 2 }($ 2 ) - (AC) 2 



7r x /4(0 2 )(d> 2 )-(AC) 5 



Kn 



vfr) 



6 2 }($ 2 ) - (AC) 2 



|AC| 



5 ( e(2Tjt2 — AC) 



(38) 



where -Ko(-z) is the modified Bessel function which can be expressed as Kq{z) = e zx (x 2 — l)~ x / 2 . This result is 
very close to the one presented in [2l| for the probability distribution function (pdf ) of the two-time composite field 
(f> a (x,i)(f>a(x,t') in the O(N) ferromagnetic model in the large N limit. In both cases, the time dependencies enter 
only through the correlation functions, {<j) 2 )(t) = C(t,t), C(t",t"), C(t',t') and C(t",t'). 



B. Composite operators 

The properties of Za™> to ) can be computed explicitly as 

n m 

V M = { Yl ^)^)JJ^.)^)/(2T)) 

i=m+l j — 1 

n rn 

= ( n mmUmafy/czT)) 

i=m-\-l j — 1 

n 

= ( II WkWk) ) £T)- m . (39) 
k=l 

The remaining average can be expanded in products of two-point correlation and linear response functions by using 
Wick's theorem applied to the Gaussian variables <fi and £. In so doing one finds the explicit 2n-time dependence of 
p(n,m)_ jf one j g i n t cres tcd in the behavior of the generalized FD ratio x( n > m ' m ' this calculation is not necessary 
since the non-trivial factors in the numerator and denominator cancel out and one simply finds a constant, 

m—m 

(40) 

independently of n. For m = 1 and rn! = one recovers X = 1/2, the usual FD ratio [HI, [UJ of the random walk. 
Equation (j40|) generalizes the result of Sec. MI Al showing that moments with different m but the same n scale in the 
same way for any choice of the time variables. 



j^(n;m,m') 



IV. THE FREE SCALAR FIELD 



The Langevin relaxation of the scalar field in d spatial dimensions is given by 

dcj)(x,t) 5F[<t>] 



m smt) +mt) - (41) 
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In the free-field case the Ginzburg-Landau functional is simply 



F[4>] = \ I 



(42) 



The expectations of the thermal noise are the usual ones reported below Eq. (p}. Equations (|4Tj) and (|42p also 
constitute the Edwards- Wilkinson model for the motion of an interface (with no overhangs) in d transverse dimensions. 
In the context of interfaces, the fluctuations of a two-time quantity the average of which is the roughness were studied 
in (23-[23. 

The Fourier transformed noise statistics are such that = and (£(£,*)£(£',*'))= 2T(2-7r) d S(k + k')5(t- 

Starting from 4>(x, 0) = 0, without loss of generality, one has 



<f>(k,t)= / dse- fc2 (*- s ^(fc,s) 



(43) 



and 4>(k,t) as well as <fi(x,t) inherit Gaussian statistics from £. From Eq. ([43 p for the Fourier space correlator 
one obtains 



(44) 



Introducing a short-distance cut-off a 2 = 1/A 2 mimicking a lattice spacing, so that J dk — > J dA; exp(— -p-), the 
real space correlation function reads 



C(r;t,t l ) = (<j>(x,t)4>(x\t')) = 
with r = |x — x'l, and one finds 



dk 
(2tt) c 



C{k;t,t') e~ lkr 



yy2-d 



1, 



4[l + A 2 (t + t')]'4[l + A 2 (t-t')] 



where 



T[n, a, 6] 



(45) 



(46) 



(47) 



is the generalized incomplete Gamma function. Analogously, the linear response is 

5{<f>(x,t)) h 



R(r;t,t') 



Sh(x',t') 



(<i>{x,m^,n)/m 



h=0 



A d 



2, .2 



4[l + A 2 (t-f)] 



(47r) d / 2 [l + A 2 (t-t')] rf/2 



9(t-t') 



(48) 



The relevant long-times limit is such that A 2 (t — t') 1. In this limit the partial derivative of Eq. (|4"6"|) with respect 
to t' becomes 



d t 'C(r;t,t') 



TA a 



l-d/2 



2 d 7T d/2 

(and similarly for dtC) while from Eq. (|48)) one has 

A d 



R(r;t,t') 



2 d 7r d/2 



t'~ d/2 (y- l)- d/2 e- c 0(i-f') 



(49) 



(50) 



where £ = r 2 /[4(t - <')] and ?/ = t/t'. Equations (gSJ and ((501) mean that for r 2 < 4(i - i'), i? and d t >C scale 
in the same way, namely R ~ t' _d/2 /ii (*/*') and <9 t <C ~ t'~ d / 2 f dc (t/t'), with / fl (y) = A d [47r(y - l)]~ d / 2 and 
/ac(y) = TA d (47r)~ d / 2 [(?/ - l)- d / 2 + (y + l)~ d/2 ], respectively. In this regime, the FD ratio (27| 



lim AT r;i,t = Hm - \' \ 
r*<(t-t') v ; r*<(t-t') d t >C(r;t,t') 



2/ + 1 



d/2 



(51) 
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is independent of r and converges, for y — > oo, to the limiting value [28| 

Xoo = lim lim X(r;t,t') = -. (52) 

Notice that this asymptotic value does not depend upon the distance for any choice of r [not only for r 2 <C (t — t')] 
since, for y 3> 1, R and dt'C are proportional in any case: 

lim X(r;t,t') = l . (53) 

t/t'-^oo 2 

This result is the same as the one found for the random walk problem, see Sec. IIIII 



A. Joint probability distribution 



In the large volume limit the joint probability (|30[) can be computed by using saddle point techniques. Changing 
the integration variables to zc = if), z^ = iX, and letting z* c = ir)*, z* = iX*, Eq. (|30f) can be cast as 



r ioa+z* » ^ ^ 

P(AC V) = I — ° / ZZE e Vh(d c ,d^,S,x' ;t,t' ,t" ,z c ,z„) 



loo+z c dz c 



ioo-\-zt 



2 71 



ioo-\-. 



2tt 



(54) 



where = = are the correlation and response densities, and 

h(dc,d li ;x,5?;t,1?,t",z c ,z li ) = -z c d c - z^d^ + G(x,x' ,t,t' ,t" , z c , z M ) 



(55) 



with 



and 



G(x,x',t,t',t",zc,z ft ) 



dk 

2 y (27r)< 



In If (fc, j?, i, i , £ , z c , z M ) 



(56) 



H(k, x,x ;t,t ,t ,z c , Z/j) = 1 - 2 Zce 



lk\x—x' 



i&ja;— a; | 



(2371) 



1 



2ife|;c— a? | 



[2T/2AC - 2C^C 9S ] + —z 2 c e Mk \ x - x [(AC) - 4C^C**] 



+ _ z 2 e 2ifc|£-S'| [(2T ~ )2 _ 4C ^ CsE 



(57) 



In the large V limit, to any choice of the fluctuating correlation and response dc and d^ there corresponds a 
couple of real quantities zq = Zq and z^ = z* whose contribution dominates the whole double integral in Eq. (|54|l . 
Zc(x,x",t,t',t") and z*(x,x',t,t',t") are solutions to the coupled system of equation: 

i — 



1 f dk ^e ik ^- s '\AC- ie 2lfe l^'l[2r^AC-2C^C <[ , H ]z;: 



i e 2^|x-x'| [(AC) 2_ 4C ^ c$$]z£ 



2 7 W H(k,a>,x>,t,t>,t",z* c ,zD 
dz/j, 

1 /• rffc i ei fc|s-£'l2T7i - \e 2 ^ s ~ s '\[2TJiAC - 2C w C7* B ]*o - | e 2i ^^-^'l [(2T/I) 2 - 4C^C ss ]z* 



(58) 



(2tt) c 



-ff (A;, x, a;', i, V , i", z^, £*) 



Once these equations are solved to find Zq and z* the joint probability distribution for large V can be written as 



1 



(59) 
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In order for Eq. (|56|) to be defined and the whole procedure to be meaningful the saddle point solutions must satisfy 
the constraint H(k, x, x'; t, t', t", z c , z*) > 0. For any choice of k, x, x'; t, t', t" , this defines the interior of the ellipses 
H = in the zj,z* plane. Since momenta are integrated over, the constraint must be obeyed for all the values of 

k. In order to sec which is the momentum which provides the most stringent condition we must know the expression 
of all the correlators entering H in Eq. ([57]) . that we derive below. Using Eq. (|43| and the properties of the thermal 
noise, with the definitions of the momentum-space correlators given below Eq. (|28[) , one readily finds 



Jl(k;t,t",t') = — 



e -k 2 (t-t") _ e -k 2 (t~t') 



C ss (k,t',t") = 2T(t"-t') 



C^ s {k,t',t") 



2T r 



1 - e 



-fc 2 (t"-t') 



(60) 
(61) 
(62) 



The other Fourier-space correlators are obtained using Eq. (j44|) and they read 



AC(k,t,t',t") = C(k,t,t") - C(k,t,t') 

k 2 i 

CW(M',0 = C(k,t",t") + C(k,t',t')-2C(k,t',t") 



e -k 2 (t+t") _ e -k 2 (t-t') 



-k 2 (t+t') 



C<f,<f,{k,t) 



T 

k 2 

C(k,t,t) = ^ 



-2k 2 t" e -2k 2 t' 2e _fe 2 (t"-t') + 2e -k 2 (t"+t') 



-2k A t 



(63) 

(64) 
(65) 



Let us notice that for k — > all these quantities converge to the value 

Ji(k;t,t",t') = C S s(k,f,t") = C^{k,t',t") = AC(k, t, t', t") = C^(k,t',t") = 2T(t" - t') = 2T5 (66) 
except Csd, that tends to the expression 



C^{k,t)=2Tt. 



(67) 



In the limit of large t — t", with the help of the expressions (|50|) - ((63)) . it is easy to check that, as t grows, the ellipse 
H = shrinks. Moreover, while for any finite k this curve approaches an asymptotic finite size as t —> oo, for k = 
the ellipses shrinks to zero. This is due to the large-i divergence of C^^, Eq. ([57]). Hence we conclude that, as k 



is varied inside the integral in Eq. (f5"7} . the most severe constraint H(k,x,x't,t',t", 



-Ci ^ 



> is provided by the 



zero momentum modes, for large t. We shortly denote with Ho(t, t' , t" ', Zq, z*) the value of the function H for k ~ 0, 
namely Ho(t, t' , t" , z^, z*) = H(k = 0,x, 5?;t,t',t", Zq,z*), and indicate with zp and the values of z* c and 



for which the constraint is satisfied H?,(t,t' ,t" , zq, z^) = 0. Let us now go back to the saddle point equations (|58[) 
Since the fluctuating quantities dc and d M appear explicitly on the left-hand-side, while z* and z* are involved into 
a complicated function on the rhs it is easier to consider the latter as independent variables, trying to find the values 
of dc and d^ for any given couple z* and z*. Approaching the constraint Hq = 0, since zc and z^ are finite and 

the denominators on the rhs vanish at k = 0, the integrands diverge. Since the first small-fc corrections to the k = 
results (|6"6")l and (|S7|) arc of order k 2 , the integral diverges for d < 2 and converges otherwise. 

In d < 2 the saddle point equations (|58|) imply that, on approaching the manifold Hq = 0, dc and d M must diverge as 



well. Reverting the argument, for dc and d^ large enough (positive or negative), the saddle point solutions z*,z* take 



nearly constant values z* 



Let us recall now that for t — > oo the size of the constraining ellipse Hq = 



vanishes, thus implying also zc — > and z M — > 0. Hence, in this large time limit the solution z c , z* to Eqs. ([55)1 is 

approaching the value zc and z^ in an ever increasing range of dc , d^ which is moving closer and closer to the average 
values dc,d^. In this range, the integrals in Eqs. (|58j) . and hence all the physics of the problem, are dominated by 



the k = behavior of the momentum space correlators, Eqs. (|66| and (|67|) . make the joint probability (|59[) symmetric 
under the exchange ~zcdc & 'z^d^, or, equivalently z^AC <H> 2Tz M /I. This is an analogue situation to the one 
encountered in the simpler case of the random walk, but now this property is only obeyed asymptotically for large 
t — t". Accordingly, one concludes that AC and 2Tfl are equally distributed in a range ever increasing with t — t", in 
any dimension d < 2 and for any choice of (x, x'; t, t' , t") provided t — t" is large. This suggests that all the composite 
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operators of Eq. (|31[) with the same value of n [with (xi,x' i ;ti,t' i ) = (x,x';t,t r ) Vi] have the same spatio-temporal 
scalings. 

Let us consider now the case d > 2. Here the integrals on the rhs of Eqs. (|58|) remain hnitc as the manifold Hq = 
is approached. This implies that, upon increasing the absolute value of ded^ up to certain finite values, a frontier F 
is met, where the limiting saddle point solutions 1c, are reached. Outside F there is no solution to Eqs. (|58]l in this 
form. This signals that the large- V limit, that we have taken from the beginning by replacing sums over momenta with 
intergrals, namely V^ 1 ~* I dk/(2ir) d , must be reconsidered more carefully. Singling out the largest contribution 
to the integrals, which comes from the k = components, in place of Eqs. (f58|) one obtains 

da = 



df, 



T8 l-± I [(8-2t)z* ll + {8-U)z* c ] 
2V H (t,t',t",z* c ,z;) 
+ 1 , y^ s ~ s '\AC - y^\s-s'\ [2TjuAC - 2CmC* s ]z* - I e 2l ^-^I[(AC) 2 - AC H C^]z* c 

2^ p H(k, x, x';t,t' ,t" , Zq, z^) 

T61-U S - 2t ) z *c + (S-M)z;} 



2V H {t,t',t",z* c ,z*) 

1 x ^ , le l ^ s - 3 '\2TJi - \e 2l ^ 3 - 3 '\[2Tji&C - 2C^C^\z* c - ^ e ^s'\ [(2TJi) 2 - 4C^,Gsb\z^ 



2V H(k,x, x';t,t' ,t" , Zq, z£) 

where the first term on the rhs is the k = term and denotes the sum over all the wavevector excluding k = 0. 
Inside F the first term is negligible and taking the large-]/ limit one recovers Eqs. (|58[) which admit a solution. 

Outside F, requiring the existence of the solution, in the large- V limit the first terms must equal dc — dc and c? M — 

respectively, while the sums transform back to the converging integrals of Eqs. (|58|) . The saddle point solution 

outside F is therefore sticked to the limiting value zc, ~z^- Interestingly enough, this implies that P(AC, ju) has a 
singular point (a discontinuity of a derivative) on F, a feature already observed in other non-equilibrium probability 
distributions [29j]. For values of ded^ well outside F, namely for large fluctuations, the contributions provided by the 
k = momentum dominate in Eqs. f|68[) . If we reason now as done for the case d < 2 we find the same conclusion as 
regards the distribution of fluctuations (namely AC and 2Tjl are equally distributed) and the scalings of the momenta. 

Let us emphasize that the scaling properties found in the large t — t" sector are fully determined by the k = slow 
momentum in any dimension d. 



B. Composite operators 

The fluctuating two-point operator the average of which is the linear response (in real space) is Ri = 
cf>(xi,ti)^(x' i ,t' i )/(2T). In consequence, the higher order correlation I?( n > m ) is given by 

n rn 

pK-) = ( ft ^,i i )0(x i) iOn^'*i)^'^)/(2T)) (68) 

i— m+1 j — 1 

This is a product of n Gaussian fields, more precisely, n — m factors (jxf) and m factors <f>£. Wick's theorem allows 
us to factor this product into products of two-field averages of the form (<fi(f>), (4>4>), (4>4>), (4>0 and (<t>0- These are 
simply C, dC, ddC, 2TR and 2TdR (where d indicates a time-derivative and one has to be careful about which is 
the time it is acting upon) . It is not difficult to sec that in the regime of largely separated times such that all ratios 
are order one 

l«A 2 (t fc -^) and ^=0(1), (69) 

ti 

and for short distances 

\x k - xi\ 2 <^ {t k - ti) , (70) 
so as to make the expressions simpler, the correlation £>(™' m ) scales as 



i=m-\-l j — 1 
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; Kn+2 \ J-d t <t™ 

-JC 7 *m+2 JddC , 

hn+1 / V *v 



fc=i 

where the proportionality is given by a function of all ratios of times. This implies that the generalized FD ratio is 
also finite 

T)(n,m) / i \ m—m' 

X)(n,m') V ' 

where the proportionality is also given here by a function of order one that depends on all ratios of the times involved 
in the 2?'s. This result is akin to the one in Eq. (|4T)|) that was obtained for the random walk. 



V. THE BIDIMENSIONAL CLOCK AND XY MODELS 



The p-state clock model is defined by the Hamiltonian 

H [a] = - J Wj = ~ J Yl cos (^ - & ) ' ( 72 ) 

(ij) (ij) 

where <ii = , <7j ) is a two-components unit vector spin pointing along one of the p directions arctan((rj /o^ ) = 
4>i = lixnijp with rii G 1,2, (ij) denotes nearest-neighbor sites i,j on a, in our case, square lattice in spatial 

dimension d = 2. This spin system is equivalent to the Ising model for p = 2 and to the XY model for p — > oo. For 
p < 4 the clock model has a critical point separating a disordered from an ordered phase at T = T\ . For p > 5 there 
exist two transition temperatures T\ and T2 > T\ [30(. For T <T\ the system is ferromagnetic, and for T > T2 it 
is in a paramagnetic phase. Between these two temperatures, for T\ < T < T2, a KT phase [311 ] exists where the 
correlation function behaves as G eq (r) ~ |r| _, '( T ) with the anomalous dimension ry(T) continuously depending on the 
temperature. Both transitions are of the KT type, namely the correlation length diverges exponentially as T\ or T2 
are approached from the ferromagnetic or the paramagnetic phase, respectively. The lower transition temperature 
goes to zero [3(J H3 (approximately as T\ ~ p~ 2 ) as p grows large, whereas T 2 remains finite. 

In the following, dynamics are introduced by randomly choosing a spin and updating it with the Metropolis 
transition rate 

w([a] -> [a']) = min[l,exp(-AS/T)] , (73) 

where [a] and [a'] are the spin configurations before and after the move, and AE = H[a'] — H [a]. In the limit p — > 00, 
in which the angle becomes a continuous variable, Langevin dynamics can also be used. We give our conventions for 
these dynamics in Sec. IV Al where we introduce the spin- wave approximation of the 2d XY model and we develop our 
analytic results. The numerical ones of Sec. IV Bl follow rule (f73|) . 

We will consider the non-equilibrium process in which a system of infinite size, initially (at t = 0) in equilibrium 
at temperature Tj, evolves for t > in contact with a thermal bath at a new t emp erature T. Various aspects of 
the kinetics of the model after such a thermal jump have been considered in [H, H3j. In the present paper, we will 
always restrict the discussion to p > 5 and final temperatures T in the KT phase, and we will consider the heating 
process starting from Tj = or the quenching protocol where Tj = 00. The heating case with p — > 00 can be treated 
analytically in the spin- wave approximation. This will be the subject of the next section. The results of this approach 
will prove to be useful also for other cases, namely quenched or heated systems with arbitrary p > 5, that will be 
studied numerically in Sec. IV Bl 



A. Heating from T = in the p — } 00 model 



In this Section we study analytically the Langevin dynamics of the clock model in the limit p —> 00, i.e. the 2d XY 
model. We first recall the known behavior of the averaged two-point and two-time correlation and linear response 
functions. Next we present our results for all moments of the fluctuating quantities the averages of which yield the 
usual correlation and linear response. 
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1. The averaged correlation and linear response 
In the spin- wave approximation the free energy functional reads (35| 

m = ^ J dS [Vmf , (74) 

where p[T) is the spin-wave stiffness. The dynamics are described by the Langevin equation 

dcj)(x,t) 5F[(f>] 



dt S<f>(x,t) 



+ Z@,t), (75) 



where the thermal noise obeys (£(£,*)) = and (£(x , t)£(x', t')) = 4Trr,(T)p(T) S(x - x')S(t- f), and the relation 
2nr](T)p(T) = T holds To ease the notation we set p{T) = 1; indeed, it is clear from Eq. (j?5j) that the actual 
behavior with p(T) ^ 1 can be recovered at the end of the calculation by a re-definition of rj(T) and a trivial time 
re-scaling. Similarly, we set fcs = 1. With these propositions Eq. ([75)) is equal to Eq. (|4ip . so that we can borrow 
the results of Sec. IIVI whenever it will be needed to infer the properties of the actual XY system. In order to avoid 
confusion between the two models, quantities relative to the scalar held will be denoted with an index ^ (e.g. and 
i?0 will be the correlation and response of the scalar held, already given in Eqs. pp]) and (|48p). 
From the knowledge of the angle dynamics the spin correlation 

C(r; t, t') = (cos[0(.f , t) - </>(£', t')]) (76) 

can be readily evaluated. The spin linear response function 



/8=1 0=1 K ' ; h=0 



(77) 



where the vector h = (h^\ h^) is the perturbation conjugated to a [i.e. by adding —p(T)/2 J dx 3{x)h{x,t) to the 
free-energy] can be obtained from 

R(r;t,t') = ±(Z(x*,t')sm[<j)(x,t) ,*')]) • (78) 



The averaged quantities C(r;t,t') and R(r;t,t') and their relation have been studied in [14l. l37H39j . We will recover 
these functions as special cases in Sec. IV A 31 



2. Composite operators 

In this paper we are interested in the more general problem of the pdf of the fluctuations of the correlation and 
linear response. In this case it is convenient to construct the pdf by evaluating all the moments. Let us start by 
defining the fluctuating quantities we are interested in as 

d = C(x i ,x' i ;ti,t' i } = cos5i(xi,^;ti,^) , (79) 
Ri = Rixi^^U,^) = — sin Si(xu 4; > ( 80 ) 

where 5i{xi, x'f, U, t'^) = (j>(x%,ti) — </>(a^,^) and, as before, x\ and x' t are two generic points in space and ij and t[ two 
generic times (still with ti > t[). In the following we will not write the explicit space and time dependences in 5- t 
to simplify the notation. For the same reason we will also set T = 1. Starting from these definitions one can build 
averages of composite fields of the form Eqs. (f3~Tj) and (|3"3"|) . These can be computed with the help of the generator 
C( n , which is obtained by adding the extra angle a\ t to <5;, namely replacing Si with = Si + a\ i in Eq. (|79p . 
More precisely, we define 

n 

4 n ' o) =<n<5i,A> (si) 
i=i 
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where 

<;.\ <■>»»,. \ (82) 



We choose the extra angle such that aA,=o = 0. Then one trivially recovers the quantity in Eq. (|33|). for the special 
case m = 0, as = d"' 0) , where {A} = means A* = Vi. Defining dC { "?l n = d Xl -d x d™' 0) 

{A}=0 W— u 

where as before d\ = d/d\ means the derivative with respect to a generic argument A, one has 



{A}=0 



dC {x}lo = ( fidaxtSmSi ) , (83) 

i=l 

where da^ = d\ i a\ i |a 4 =o- This quantity provides all the composite fields in Eq. pip when the choices ot ® for 
i = 1, ... ,771, and ck for i = m + 1, . . . ,n, are respectively made, with 

da { ^ ] = -5^(4,4), (84) 

flag) = • (85) 

The computation of the generator (|8~Tj) (see App. yields 

C (n,0) =2 -„ £ e -^(5S> 5 (86) 
{«i=±l} 

where 5a = S"=i s i&i,\ and s j are auxiliary Ising variables introduced for convenience. In the long times sector in 
whi ch 2A 2 (t + t') > 1 and 2A 2 (i - i') > 1, so that we can use the limiting behavior of and discussed in 
Sec. IV A 11 one obtains 

9C^°2 = 2-« £ (n^pW^^De^^, (87) 

{s s =±l} \*=1 / 

where 5 = 5a=o and the polynomial is defined by the recursive relation given in Eq. (|A3|) . Given the form of S 
and Eqs. (|84[) and ([55]) these quantities arc all expressed in terms of the angle correlation and responses, and Rs, 
and can, therefore, be explicitly calculated. 

Let us come now to the composite operators 2?( n > m ) defined in Eq. ([3"Tj). which can all be obtained from Eq. ((57)) . By 
comparing two different moments with the same n but different m (say m and m! > m) the difference is only due to 
the replacement dce^ — > daf C ^ Vi = m+1, . . . , m' in the polynomial ?W, which in turn amounts to the substitution 
of the functions dfC^Xi — a^|; with R^Qxi — t i: t ■) Vi = m + 1, . . . , m' . For \x% — <€. 2A 2 (i, — i-) Vi these 
quantities arc proportional, according to Eq. (|51[) . Therefore one concludes that all the composite operators (and 
hence all the moments) with equal n scale in the same way, and the generalized FD ratios (|3"4")) depend only on the 
ratios = ti/tt and are independent of the distances |af$ — x'^. Furthermore, according to Eq. (f52|) . for j/j — 5- oo Vi 
the limiting value 

X (n; m ,m') _ x ( n;ro ,m') (88) 

t/t'— >oo 

is finite and independent of the spatial arguments. 



Scaling of composite operators and their FD ratio for n = 1, 2 

As a concrete example we now cxplicitate the expressions obtained in Sec. IV A~2l for the simplest cases with n — 1, 2. 
The generalization to generic values of n is straightforward. Starting with the case n = 1, from Eq. (|86j) one 
immediately obtains 

£(i,o) = c(r;t,t') = e - '^'*)-"^''*')] 2 = e -£[c<»(o;t,t)+cv(0;t',t')-2C^(r;t,t')] (gg) 
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where we dropped the sub-index 1, namely we set Xi = x and similarly for the other variables. For the computation 
of the 2?'s we enforce Eq. (|57|) using, according to Eq. (|A3[) . 



pW({(Sda Xi }}) = -s 1 (S 1 da Xl ) 



Hence one arrives at 



and from here 



= d t ,C(r; t, t') = [d t ,C^(r- 1, t') - 9 t ^(0; i', i')] C(r; t, t') 



and 



D^ 1 ) = iZ(r; t, t') = [Rt(r; t, t') - i? (O; t', t')] C{r; t, t') . 

Assuming R^(r;t,t) = dC,j,(r;t,t) = 0, for the FD ratio X(r;t,t') = X^ 1 '^ (r;t,f) one finds 

X{r;t,t') = X^tJ) , 

showing that the FD ratio of the XY model is the same as that of the scalar field. 

Let us now consider the case with n = 2. Proceeding analogously to the case n = 1, from Eq. (|86|) one has 



c (2,o) 



: 



(90) 
(91) 
(92) 
(93) 
(94) 

(95) 



where 



4 2 ' 0) = cxp|-i[C (O;i 1 ,t 1 ) + Q(O;t 1 ,t 1 ) + C (O;t 2 ,t 2 ) + C4O;4,4) 

-2C4\x 1 -x[\;t 1 ^[)-2C^\x 2 -x' 2 \;t 2 ,t' 2 )}} 
x exp{±2[C !j) (\x 1 - x 2 \;t 1 ,t 2 ) - C 4> {\x 1 - x' 2 \;t 1 ,t' 2 ) 

-C (|fi - x 2 \;t 2 ,t' 1 ) + C^(\x' 1 -x 2 \; t[ , t' 2 )\ } . 

The second order recursive polynomial reads 

P^{{{Sda Xi )}) = (Sda Xl )(Sda Xl ) - s 2 (da Xl da X2 ) 



(96) 



(97) 



where, using the terminology of App. [A] the two terms on the rhs are of type / and /, respectively. Neglecting the 
latter, since we show in App. [B]that it is sub-dominant in the large time sector, we arrive at 



a r (2,o) l 
oc {A}=o - 2 



((<*! + S 2 )da Xl )((S 1 + 5 2 )da X2 )cf a) 
+<(*! - <S 2 )da Al )((<$! - S 2 )da X2 )C^ 0} 



or, more explicitly, 



with 



2 



A (2, m ) c (2fi) + A (2,m) c (2,0) 



(98) 



(99) 



.4 



(2,0) 



[d t >C <l> (\x 1 -x' 1 \;t 1 ,t' 1 ) - dfC^t'^) 

±d t > i C <l> (\x 2 -xi|;i 2 ,ti) Tdt^C^x^ - x' 1 \;t' 2 ,t' 1 )] 

±dt'C <p (\x2 - x' 2 \;t 2 , t' 2 ) T drC^O; t' 2 ,t' 2 )} , 

4 2 ' 1} = [d t ,C^\x 1 -x\[MA)-dt\C^t' 1 A) 

±d t > i C <j> (\x 2 - x* x \;t 2 ,t' x ) T d^C^dx^ - x' 1 \;t' 2 ,t' 1 )] 
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x [R^\x 1 -^\-,t 1 ,t' 2 )-R (j> (\x' 1 -x' 2 \-,t , 1 ,t , 2 ) 

±R (j> ( \x a - 4 1 ; t 2 , f 2 ) T #0 (0; t' 2 , t' 2 )] , 
A [ l> 2) = [R 4) {\S 1 -x' 1 \;t 1 ,t' 1 )-R^{0;t' 1 ,t' 1 ) 

±R 4> (\x 2 -x'^-hA) T^(|4 
x - 4|;ti,t' 2 ) - - 4l;*i>*2) 

±i?. (|f 2 -4|;t 2 ,4)T^(O;t' 2 ,t' 2 )]. (100) 

In the following, in order to simplify the discussion, we focus on the case x\ = x\ ,x' 2 = x 2 and ii = ti = 4 , t'-y = 
t' 2 = t'. This choice will be adopted in Sec. IV Bl for the numerical computations. Letting r 2 = (xi — x 2 ) 2 <C 2A 2 (t — t') 
and using the scaling form for R^ and dfC^ derived in Sec. QV] [below Eq. ([50))]. one easily obtains 

A ;( n =2;m,m')( r;fjt /) = [^(y)]™"" 1 ' , (101) 

where y — t/t' as usual. The limiting FD ratio is given by 

( 1 \ rn — m' 
• (102) 
2 J 

Proceeding analogously, one can generalize the computation of X^ m ' m ' to any value of n, m, ml (but the calculation 
becomes lengthy upon increasing n). 



4- Moments 



Moments can be defined from the composite operators (|31j) and ([33)1 by subtracting a suitable disconnected part. In 
view of the numerical applications of Sec. IVBl we will concentrate here on the time- integrated quantities of Eq. (|33l) . 
These quantities are less numerically demanding than the corresponding differential ones pip . For the same reason 
we restrict to the case n = 2. With the choice x\ = x\ , x! 2 = x 2 and ti = t 2 = t ,t[ = t' 2 = t' made in the previous 
section, moments can be defined as 

V {2 ' m+m '\x u x 2 -tA) = (^• m \x 1 ;t,t')C^ 1 ' m 'Xx2;t,t')} -C^ m \t,t')C^ m 'Xt,t') (103) 

where C is the fluctuating part (insides brakcts (. . .)) of C in Eq. (|3"3"|) . In order to improve the statistics of the concrete 
numerical measurements that will be discussed in the next section, and to compare to similar calculations presented 
in (ill . Il2j , we compute the double spatially integrated quantities 

V&frO = L ~ 2 I d ^ I d ^ V^ m \x x ,x 2 ;tA) (104) 



with L the linear size of the sample. y( 2 >°) is the quantity that is usually computed when dynamical heterogeneities 
in disordered and glassy systems arc studied [4(J ■ Together with this one, the other two quantities have been studied 
in different aging systems and with different techniques in 

Using Eqs. (jHFJ and V^Jq 1 can be written as 

y(2.0)^ _ l e -C^(0;t,t)-C^(0;t',t')+2C^(0;t,t') 

dr r e 2 C*(nM)+2C^(r;i',i')-4C4,(r;t,i') e -2C^(r;t,t)-2C^(r;t' ,t')+iC^(r;t,t') 

where r — x! — x. Using the expressions derived for in App. [C]for 2A 2 (t - f') > 1 wc obtain the scaling form 

v£2 ) =* ,0 / {3 '%) (105) 

with a = (2 — 2rj)/z, where r\ = T/(2n) is the equilibrium anomalous exponent and z = 2 is the dynamical exponent. 
This result agrees with the general behavior 

v (2,m) = t , { A- d -2 n) / Zf (2, m ) {y) (1Q6) 

expected on the basis of critical scaling arguments 11 1, and confirmed in (loj in the spherical model. The same scaling 



is obeyed also by V^_Lq and VL^q , with the same large-i/i' behavior of the scaling functions, since we have proven 
in Sec. IV A 21 that all composite operators with the same n scale in the same way. Notice that the exponent a is an 
equilibrium property, being only determined by the equilibrium exponents ry and z. 



19 



B. Numerical simulations 

The analytical results of the previous section apply to a system with p — > oo heated from T, = to a temperature T 
in the KT phase. The next question is how general this picture is and, in particular, i) what is the behavior of systems 
evolving in a KT phase with p < oo and ii) which modifications arise in a quench with Tj = oo where topological 
defects are present due to the disordered initial condition. In this section we address these questions numerically. 
In order to do so we evolved systems with p = 6 and 12 starting from equilibrium states at Tj = and Tj = oo. 
In the case of a quench in the p = 6 case, the behavior of C(0;t, t') and x(0;4, t') was shown (34} to fit into the 
general scenario expected from standard scaling arguments (apart from logarithmic corrections due to the presence 
of vortices, see the discussion below). Here we will concentrate on the behavior of the moments of Eq. (|104[) 

which, in the present on-lattice model are obtained as V^f Q (t,f) = L^ 2 ^2 i ^-V^ 2 ' m \xi,Xj]t,t l ), where Xi (xj) is 
the (square) lattice coordinate of site i (j), and L 2 is the number of lattice points. Whenever a response function 
is involved this has been computed with the extension of the FD theorem to non-equilibrium states derived in [4lj . 
This method has been thoroughly applied [42[ to study different problems for its numerical efficiency and because, 
being perturbation-free, guarantees correct results in the linear regime. The working temperature is chosen to be 
T = 0.76, which belongs to the KT phase both for p = 6 and p = 12, and the system size is L = 600. No finite size 
effects arc detected with this choice, in the range of simulated times. The data presented are averages over 2 • 10 3 - 
6 • 10 3 (according to the different cases) realizations of the thermal noise and, in the case of quenches, of the initial 
conditions. Times are given in Monte Carlo (MC) units. Whenever we plot the moments, we include a suitable T 

factor to make them dimensionless; namely we always plot v£^q\ TV^q , an d T 2 vj? = l^ > . 

1 . Heating from zero temperature 

Clock model with p = 6. 

The numerical estimates rj = 0.17 and z = 2.18 are reported in the literature [32l . l34j |. from which, comparing 
with Eqs. ([105]) and (fT06|) one obtains a = (4 - d - 2r])/z ~ 0.76. In Fig. Q] the quantities t~ a V^L^ are plotted 
against t/t' in the upper panel. One observes a nice data collapse for V^Iq and V^J , where small corrections are 

(2 2) 

only visible for the smallest value of t in the early regime with t/t small. The quantity V^Lq, instead, presents 
larger corrections, and only an asymptotic trend towards a scaling collapse is observed (the two largest value of t' 

(2, 2) 

are almost superimposed). The same large corrections to scaling affect also the scaling function of Vj._l . Indeed, 
while f( 2 '°>(x) and f^ 2,x '(x) are proportional, as expected, and grow algebraically as f^ 2,0 \x) ~ Z^ 2 ' 1 -' ~ x a , with 
a value of a ~ 0.8 (measured for x > 10), f( 2 ' 2 \x) grows with a larger effective exponent ot e ff but, as it is more 
visible for the largest t', this exponent decreases as x increases (for t' = 100 one measures a e ff — 0.9 for x > 70). 
In order to better appreciate the fact that all moments scale with the same exponent, and to test this fact in a 
parameter-free plot, in the lower plot of Fig. [1] we present the parametric plots of —Vj^Q and Vj^q against V^q . 
We find an excellent data collapse for the different values of t', confirming once again that all the moments scale in 
the same way. Notice also that data follow a linear behavior (green line) for large times, signaling that also the the 

(2 2) 

scaling functions are proportional (for V^_! , due to the above-mentioned pre-asymptotic corrections, the approach to 
a linear behavior is seen only for the latest data). In conclusion, our data are consistent with an asymptotic scaling 
"^fcio = £'(4-d-2r;)/z j(2,m) (t/t'), with the scaling functions increasing algebraically with an m-independent exponent 
a. The limiting FD ratios limt/^oo lim t _j. 00 X^ 2 ' m ' m \t,t') are therefore finite. We conclude that the scaling scenario 
given in Eqs. (|105[) and (|106p and suggested by the spin- wave results applies to this case, provided that the actual 
values of rj and z are taken into account. 



Clock model with p = 12 

The results for p = 12 are presented in Fig. [5J One observes the same qualitative behavior as in the case with 
p = 6. With a value a = 0.83 one obtains an excellent data collapse for and V^J , while for ^ ne collapse 

is only asymptotically approached, similarly (but the collapse is somewhat better) to the case with p = 6. Notice 
that this value of a is basically the same as the one obtained for p — > oo, where one has a ~ 0.83 since 77 ~ 0.165 and 
z = 2 [43| . From p = 12 onward, therefore, one does not expect to see any significant difference with the spin-wave 
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FIG. 1: (Color online.) Clock model with p = 6 heated from T = to T = 0.76. Upper panel: The re-scaled moment 
t'~ a Vj^ j (t, t'), with a = 0.76, is plotted against t/t' in the main part of the figure, while the same plot for V&J (with a minus 
sign in order to have a positive quantity) and V^J is presented in the upper-left and lower-right insets, respectively. Different 
choices of t' correspond to different curves, see the key. The straight green bold line is the power law (t/t') ' 8 . Lower panel: 
Same data as in the upper panel but plotted in the parametric form — V&q' (t, t') against V^J^ (t, t') (main part) and (t, t') 
against (t,t') (inset). The straight green bold line is the linear behavior (i.e. scaling functions are proportional). 



analytic results. The scaling functions behave as /^ 2,0 '(x) ~ /' 2,1 - ) <~ x a , with a = 0.84. For f^ 2,2 ^ the data are 
consistent with an asymptotic convergence towards the same power-law behavior. This picture is confirmed by the 

( 2 1^ f 2 2 ^ ('20') j_ 

parametric plots of — V^_! and V^_! against V^ presented in the lower panel of Fig. [5J 



2. Quench from infinite temperature 

When quenching from Tj — > oo the presence of vortics makes the dynamics quite different from the one observed 
in the heating case [44j . It is therefore interesting to see how the scenario provided by the spin- wave approximation 
may or may not be modified in this case. In the following we present the results of simulations of the same system 
considered insofar, but with initial conditions extracted from an Ti — > oo equilibrium ensemble. We only discuss the 
case with p = 6 since we found very similar results for p = 12. 

f2 0) f2 1^ (2 2l 

The second moments V ( k= l ' (t,f), V x k= l ' (t,f) and V k Zl'{t,t') are plotted in Fig. O In this fi gure we used the same 
scaling procedure for the heated system, with the same exponent a since, being an equilibrium quantity it should not 
depend upon the non-equilibrium protocol. As it is seen, this produces a quality of data collapse comparable to the 
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V (2 ' 0) (t,t') 



FIG. 2: (Color online.) Clock model with p = 12 heated from T = to T = 0.76. Upper panel: The re-scaled moment 
t'~ a Vj^Q (t,t'), with a — 0.83, is plotted against t/t' in the main part of the figure, while the same plot for (with a 

minus sign in order to have a positive quantity) and is presented in the upper-left and lower-right insets, respectively. 

Different choices of t' correspond to different curves, see the key. Lower panel: Same data as in the upper panel but plotted 
in the parametric form —V^2^(t,t') against V^f^ {t,t') (main part) and V&q (t, t') against Vj^l^(t, t') (inset). The straight 
green bold line is the linear behavior (i.e. scaling functions are proportional). 



case of the heated system. Strictly speaking, one should expect logarithmic correction to this scaling behavior, caused 
by the presence of vortices [44| . However, such corrections are very tiny in the asymptotic time domain and they 
cannot be detected from the inspection of our data. A major difference with respect to the heated system is, instead, 
the behavior of the scaling functions. Indeed, while j^ 2 ' 2 -* keeps growing with the same power law of the heated 
system (the two quantities - in the heated and quenched case - are basically indistinguishable except for small x), 
/( 2 '°) and J*- 2 ' 1 -* are changed to a much slower, logarithmic growth. From the parametric plots presented in the lower 
panel of Fig. [T]one argues that Z^ 2 ' -' and J*- 2,1 - 1 are still asymptotically proportional. Interesting, in an intermediate 
time regime also ,p 2 ' 2 ) grows proportionally to the other scaling functions. For larger times, however, there is a 
crossover to a faster growth. This interesting feature shows that V^ (t, t') does not feel the presence of vortices, 
while the other moments do. Since these quantities have been proposed |g, 13, [111, H(| as efficient tools to detect 
dynamical heterogeneities and to quantify cooperative lengths, our results suggest that different lengths are encoded 
in the different Vs. A possible explanation could be that vj^f^ and detect the distance between vortices while 

(2 2) 

1^=0 is only determined by the typical length of the smooth spin rotations (spin waves), but this subject should be 
further investigated. Notice also that the peculiar scaling of the Vfc=o's found in the case of the quench has non-trivial 
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consequences on the limiting behavior of the generalized FD ratios of Eq. (f34|) . Indeed, since the ordinary FD ratio X 
is finite [45[ , the disconnected terms subtracted off in Eq. (jl03[) have the same scaling properties independently on m 
and m'. Hence the scaling of the V's directly inform us on the behavior of the composite operators C( n,m ' (|33|) and, 
in turn, of the ratios X^ 71 '" 1 '™ 1 > of Eq. (fM|). Since the V's are found to scale with the same exponent a but with a 
different form of the scaling function this implies that the limiting value lim t <_ i . 0O x^ 2 ' ,m ' m \t = yt',t') (with y fixed) 
is finite. The same is true also when the limit lim t _ i . 0O X^ 2 ' m ' m ^(t,t') is taken (with tl sufficiently large) but only for 
ra^2 and m! ^ 2. On the other hand, the same quantity is not finite for m = 2 or m! = 2. This behavior is radically 
different from all the other cases studied insofar. 





V (2 ' 0) (t,t') 



FIG. 3: (Color online.) Upper panel: Clock model with p — 6 quenched from T = oo to T = 0.76. The re-scaled moment 
t'~ a Vj^Q (t, t'), with a = 0.76, is plotted against t/t' in the main part of the picture, while the same plot for V^Lq (with a 
minus sign in order to have a positive quantity) and Vj^Q is presented in the upper and lower insets, respectively. Different 
choices of t' correspond to different curves, see key in the figure. Lower panel: Same data as in the upper panel but plotted 
in the parametric form — V^Lq it, t') against V^Jq {t,t') (main part) and V^Jg (t,t') against Vffi? (t, t') (inset). The straight 
green bold line is the linear behavior (i.e. scaling functions are proportional). 
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VI. CONCLUSIONS 



In this paper we undertook the study of the out of equilibrium dynamics of some unfrustrated models with a finite 
FD ratio from the novel perspective of dynamic fluctuations. After defining in a proper way the fluctuating quantities 
the average of which yield the usual two-time correlation and linear response function, we evaluated the properties 
of their probability distribution and the scaling of the composite operators [defined in Eqs. (|31j) and ((33)) ] made of 
products of n of these fluctuating quantities. We showed that, in the model studied analytically, such composite 
fields in the asymptotic time domain scale in the same way when the total number n of fluctuating parts involved 
is the same, irrespectively of how many factors are of the correlation or the linear response type. Therefore ratios 
between composite operators with the same n converge to a finite value in the large time limit for any value of n. 
Since such composite operators are strictly related to higher order correlation and response functions, these ratios 
can be regarded as a generalization of the usual FD ratio above linear order. In the restricted context of the simple 
models considered analytically in this Article, their finite asymptotic value might speak about the significance of the 
notion of an effective temperature associated to the FD ratio. Indeed, for such a concept to be physically meaningful, 
one would ask such effective temperature to remain finite at any order. Related to that, the analytical results of 
this Article support the idea that fluctuations in aging systems are intimately related to the behavior of the FD 
ratio 0, Q ■ The mechanism whereby a finite value of the generalized FD ratios is attained in the simple unfrustrated 
models considered analytically here might also be useful to understand the behavior of fluctuations in more realistic 
systems. 

Besides, the analysis of this Article is also related to the problem of the detection of characteristic lengths form 
higher order correlations and response functions 0, [^, One major problem in this context is which lengths arc 
these higher order quantities sensitive to, and how. Interestingly enough, in all the models considered in this Article 
where a unique growing length is present the composite operators have the same asymptotic time-scaling (exponents 
and scaling functions) . This indicates that the growing length enters different composite operators in the same scaling 
way. The different scaling functions found in the quenched clock model, where two different lengths are present, seems 
to indicate that different composite operators are sensitive to different lengths. 

Some previous studies of the second-order momenta in different model systems had been performed and we wish 
now to confront our findings to these. In so doing, we will refine the picture of the dynamic scaling of fluctuations 
in i) coarsening systems quenched below their critical temperature; ii) disordered spin models such as the Edwards- 
Anderson spin-glass; iii) critical systems as the ones we studied here. 

The results for the critical cases considered here (except possibly from the quenched clock model) are clearly different 



from what has been found in sub-critical quenches of simple coarsening systems [11| where the moments associated 
to the correlation scale differently from the ones where the fluctuations associated to the linear response enter. The 
analysis of the second-order momenta for a ferromagnet quenched to its critical point computed numerically in [TTJ] 
and analytically in the spherical model in [Toj unveil a behavior analogous to the one found in this paper and the 
scaling found by these authors conforms to the form given in Eq. (|I06p . Monte Carlo simulations of the sub-critical 
dynamics of the 3d Edwards- Anderson spin-glass [f2[ suggest that the second moments in these glassy systems scale 
in the same way, in agreement with the conclusions arrived at in [46l - l49j by analyzing the joint probability distribution 
of Sec. Ill Al This claim has to be taken with the usual proviso that numerical simulations of glassy systems are hard 
to interpret beyond any doubt. 

The analysis in this paper could be applied to other systems with glassy dynamics; thus helping to complete a 
general comprehension of fluctuations in problem with slow dynamics. Obvious candidates are kinetically constrained 
models [5l| , for which an analysis of out of equilibrium fluctuations along these lines was initiated in (52j ; the one 
dimensional Glauber Ising model [53j . or random manifold problems, studied from an averaged perspective in [lij 
among many other papers. 

Appendix A: Calculation of C^ 1 ' ^ 

Using the exponential form of the cosine we have 

n 

4"' 0) =2-«(JJ(e^+e-^))=2-" ]T (e*^.- Sl5 " A ) , (Al) 
»=1 {s;=±l} 
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where Sj = ±1 arc sign variables. Using the property (exp(±ia)} = cxp(— (a 2 }/2), holding for any linear function a(£) 
of £, since the Si^\ are themselves linear, one arrives at 

C (n,0) =2 -» £ e -K(E i=1 ,„^) 2 > =2 -n £ e-^), (A2) 
R=±l} {si=±l} 

where the last equality defines the function S\. We are now able to compute the n-times derivatives involved in 
Eq. ([3T]) . By taking them one at a time it is easy to check that d\ 1 ...d\ n e~ 2 < 5 a> — (J|" =1 Sj) p( n ) (5a, {daA; })e~2 
where the polynomial P" can be obtained by the recursive relation 

P^(S x ,{da Xi }) = -P^iSxAdax.mSxdax^ + s^P "" 1 ^, {dc^}) , (A3) 

starting from P(°'(S\, {da Xi }) = 1. Hence one has 

*#£>o = 2- £ fn^PW^^We-*^), (A4) 

{s 4 =±l} \i=l / 

where 5 = iS{a=o} ■ The recursive equation (|A3| implies that P^™' is a sum which contains only products of correlators 
of the type 

f=(Sda Xi ), (A5) 

and 

/ = (da Xi da Xj ) . (A6) 

Indeed, the first term in the rhs of the recursive equation (|A3[) is itself a multiplication by a factor of type /, while 
the second term amounts to the replacement of some terms of type / with others of type / (because sj d Xj (Sda Xi ) = 
(da Xi da Xj )). It is easy to show (see App. [Bj that in a large time limit [for t' — > 00, with y and r (or £) finite, or 
equivalently for t — t' — > 00 with r (or £) finite] the terms of type / are sub-dominant. Hence, in this time sector one 
has PW ~P^({(Sda Xt )}). 



Appendix B: Asymptotic analysis 

Terms of kind / are sums (over i) of contributions of the form 

f{x i ,x' j ;ti,t' j ) = (0(xi,ii)9 t /0(Xj-,^-)) = dt'.CjyQxi - x'^U,^) , (Bl) 

for j < m, or 

f{xi , x' ■ U ,t' j ) = — (tf>(xi , ti)£(a:$ , t$)> = M I ^ " ^ 1 1 *i » *J) ( B2 ) 

for j > to. Next to these there are analogous contributions that are obtained from the terms (|B1[) and (|B2[) with the 
replacement —> (x^t^). 

As it can be seen from the recursive equation (|A3[) . at the r-th step, terms of type / are generated from the 
quantities of type / already present at step r — 1 by replacing S with a da Xr . Then, in a generic term / it is always 
i ^ j. Using Eq. ([85]) for the quantities of type J with i,j > to, one has J(x' i ,x' ] \t' i ,t' j ) = (2T)~ 2 (^(x r i ,t^(x r j ,t r j )), 
which, restricting to the case, (x'^t^) ^ (x'j,tj) Vij, vanishes identically. On the other hand, for i,j < to, using 
Eq. dS3J) one has 

fix'uXjlt'iitj) = (di' i 0(^,*O% < H a; i> i j-)) = dt'dt'C^rij;^,^) , (B3) 

and for i < m and j > m 

Jix^x'^t'^) = ^(dt^t^'pt'j)) = dt.R^n^t') . (B4) 

Comparing with the corresponding terms (|B1[) and ()B2|) of type /, there is an extra time derivative in the quantities 
/. Using the scaling properties (|4T)|) and ([50)1 one concludes that the latter are negligible in the large times domain 
[for i! — y 00, with y and r (or £) finite, or equivalently for t — f' — > 00 with r (or £) finite]. 
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Appendix C: Scaling of C$(r;t,t') 
From Eq. (g6j), for 2A 2 (i - 1') > 1 one can write 

T e~ z 

C (r;t,t') = — / (CI) 

where £ = r 2 /[2(i — i')] and y = t/t'. For £ <C 1, neglecting the exponential factor one obtains 

C 4 (r;t,t') = ^]n(j±l), (C2) 

For larger values of £, instead, letting ?/ 3> 1 one finds 

C (r;t,O = ^(y-l)- 1 e^. (C3) 

Z7T 

Finally, we consider the equal times correlation. Proceeding as before we write 

C^r ]t ,t) = — / dz— . (C4) 

ATT J Agrg z 
4A y t+l 

For A 2 r 2 < 1 + 4A 2 i, again neglecting the exponential, one has 

CV(r; i, i) = ^lln(l+4A 2 t) (C5) 

47T 

For 1 < A 2 r 2 < 1 + 4A 2 i one can split the integral as 

^"-cl/L. ? + (C6) 



T , / 1 + 4A 2 t 



and hence 

C,(r;M)c-ln^-^j , (C7) 

since the second integral in Eq. (|C6[) is exponentially suppressed. Finally, for A 2 r 2 3> 1 + 4A 2 £, C^,(r;t,t) is exponen- 
tially small. Putting everything together one has 

T /1 + 4A 2 A / A 2 r 2 \ 

C *^^ ^(aWi) ^(tTuh) ' < C8 > 

where f cu t(z) is a cut-off function. 
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